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Large scale circulations and energy transport in contact 
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1 INTRODUCTION 



Contact binaries have been defined bv iKuipeil (|l94ll ') as bi- 
naries with components surrounded by a common envelope. 
Within the Roche approximation for the total potential in 
a rotating frame of reference, the envelope lies between the 
inner and outer critical equipotential surfaces, defined re- 
spectively by the Lagrangian points Li and L2 . Cool contact 
binaries with spectral types late r than FO are called W UMa 

type binaries (|Mochnackil 1 198 j ) . 

W UMa type stars are fairly common in space. iRucinskH 
l|2002l ) estimates that one such binary occurs per 500 main 
sequence stars in the solar vicinity but these stars do not 
appear among members of young and intermediate age clus- 
ters. Their occurrenc e rapidly i ncreases in old open and 
globu lar clusters (Kaluznv fc Rucinski. 19931 : iRucinskill 199^ . 
I2OO0I ). This result indicates that cool contact configurations 
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ABSTRACT 

A hydrodynamic model for the energy transport between the components of a con- 
tact binary is presented. Energy is transported by a large-scale, steady circulation 
carrying high entropy matter from the primary to secondary component. The circu- 
lation is driven by the baroclinic structure of the common envelope, which is a direct 
consequence of the nonuniform heating at the inner critical Roche lobes due to un- 
equal emergent energy fluxes of the components. The mass stream flowing around 
the secondary is bound to the equatorial region by the Coriolis force and its width is 
determined primarily by the flow velocity. Its bottom is separated from the underly- 
ing secondary's convection zone by a radiative transition layer acting as an insulator. 
For a typically observed degree of contact the heat capacity of the stream matter is 
much larger than radiative losses during its flow around the secondary. As a result, its 
effective temperature and entropy decrease very little before it returns to the primary. 
The existence of the stream changes insignificantly specific entropies of both convec- 
tive envelopes and sizes of the components. Substantial oversize of the secondaries, 
required by the Roche geometry, cannot be explained in this way. The situation can, 
however, be explained by assuming that the primary is a main sequence star whereas 
the secondary is in an advanced evolutionary stage with hydrogen depleted in its core. 
Such a configuration is reached past mass transfer with mass ratio reversal. Good 
agreement with observations is demonstrated by model calculations applied to actual 
W UMa-type binaries. In particular, a presence of the equatorial bulge moving with 
a relative velocity of 10-30 kms^' around both components of AW UMa is accounted 
for. 

Key Virords: stars: contact - stars: eclipsing - stars: binary - stars: evolution 



are not formed at, or near the zero age main sequence 
(ZAMS) but at the substantially later age. As it is now 
commonly assumed, they are formed from initially detached 
binaries which lose orbital angular momentum (AM) at such 
a rate that a contact is reached after several Gyr. Kinemat- 
ically, field W UMa s tars belong to old disk, which supports 
their advanced age (|Guinan fc Bratstree^ll988l : iBihr et all 
12005!). The most promising mechanism for AM loss is related 
to the chromospheric-coronal activity of bin ary components 
l|Vilhulll982l : lMochnackilll985l : IStepienlll995l ) although a high 
incid ence of companions to W UMa stars (|Rucinski et al.l 
l2007t ) suggests t hat a third body may als o play a role in 
orbit tightening l|Eggleton fc Kiselevalf200ll ). 

lKuipeJ(| 19411 ) noted that a contact binary is stable only 
when both components are identical. Otherwise, a binary 
consisting of two stars with the same composition is unsta- 
ble and mass will be transferred inside the common envelope 
from the more massive (primary) to the less massive (sec- 
ondary) component. The transfer is driven by a nonuniform 
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heating of the base of the common envelope due to unequal 
emergent luminosities at the inner critical surface. By conti- 
nuity, temperature differences on each equipotential surface 
above the Roche lobes exist. As a result, the common en- 
velope must be treated as baro clinic rather than barotropic 
l|Shu et al.lll"979l : lTassoullll992l ). Different vertical (i.e. per- 
pendicular to the local equipotential surface) pressure strat- 
ifications in both components produce the horizontal pres- 
sure gradient driving mass motions on a dynamical time 
scale. Even if the pressure is constant over a given equipo- 
tential surface, a diff erence appears above and below that 
surface (Kah terl 120041 '). In the absence of the Coriolis force 
the resulting large scale fiows are symmetric a round the axis 
joining centers of both stars (|Webbinklll973 ). However, for 
a flow velocity being a significant fraction of the sound ve- 
locity the Coriolis force cannot be neglected. The force acts 
outwards in the equatorial plane, defiects the fiow towards 
one side of the neck between the components and makes it 
go around the other star in the direction of its orbital mo- 
tion. Figure 1 shows geometry of both critical Roche lobes 
and of the fiow in the equatorial plane (see also Fig. 4 in 
iKuipeij [l9411. Subsequent computations of the streamlines 
in semidetached binarie s confirmed a strong i nfiuence of the 
Coriolis force (see e. g. iLubow fc Shul Il975l , iLubow fc Sh^ 
ll976l . lOka et al.ll2002l ). When the flow returns to the parent 
star it is directed to the other side of the neck (Fig. 1). 

Together with mass, thermal energy is carried to the 
other component, so we should expect a signiflcant mod- 
iflcation of the apparent surface brightness of both stars, 
compared to a detached binary. Indeed, observations of W 
UMa type stars show eclipse minima of nearly equal depth 
indicating almost the same surface brightness of both com- 
ponents. It indicates that, unless mass ratio is close to one, 
a substantial fraction of the flux radiated by the secondary 
component comes from the primary. 

A very elegant model of a W UMa type bi nary consist- 
ing o f two main sequence s tars was d eveloped by ^ Lucy* ("196^, 
1 19761 ) and lFlanner"^ (| 19761 ) (see also I Yakut fcE gglcton 2005 
and references therein). The model assumes that each com- 
ponent is out of thermal equilibrium, with a radius oscil- 
lating around its critical Roche lobe (inner critical surface). 
The model requires that matter transported from the pri- 
mary fully covers the secondary like a hot blanket. The blan- 
ket blocks completely the radiation energy flux produced in 
the core of the secondary. The blocked flux is converted into 
the thermal energy of the secondary which expands on the 
stellar thermal time scale. Specific entropy in its convection 
envelope increases, approaching the value corresponding to 
the primary's convection zone, which, at the same time, de- 
creases somewhat due to energy transfer. When the expand- 
ing secondary overfiows its Roche lobe, mass and energy 
transfer from the primary is stopped and the star shrinks 
radiating away the excess thermal energy until the cycle 
of these thermal relaxation oscillations (TRO) repeats. The 
TRO model successfully explained two important properties 
of W UMa type stars: an abnormal radius ratio of the com- 
ponents (assumed to be MS stars), required by the Roche 
geometry, and essentials of the observed light curves. 

In this paper, calculations are presented that reveal dif- 
ficulties with the TRO model. In particular, it will be shown 
below that the Coriolis force is strong enough to balance the 
meridional pressure gradient of the circulation current. As 
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Figure 1. Geometry of a contact binary in an equatorial plane. 
The binary rotates anticlockwise around the axis crossing the 
plane in a point marked with a cross. Broken lines denote inner 
and outer critical surfaces. A common envelope lies between the 
inner critical surface and the stellar surface plotted with a solid 
line. Mass flows from the primary through a neck (the stellar 
surface narrowing) and encircles the secondary, as indicated by 
arrows. (For a side view of a contact binary see Fig. 5). 



a result, the stream of matter from the primary is bound 
to the equatorial belt of the secondary. The polar regions 
are not covered by the stream so the star can freely radi- 
ate away its nuclear energy. As a result, the specific entropy 
in its convection zone increases insignificantly and remains 
considerably below the value characteristic of the primary 
convection zone. Also its radius hardly increases, which ex- 
cludes MS secondaries oversized by a f actor of 2-3 observed 
in many W UMa stars (|Stepiehl l2006al ) . Substantially over- 
sized secondaries can, however, be naturally explained by 
assuming that, instead of being MS objects, they are highly 
evolved stars with depleted h ydrogen in the center o r even 
possessing small helium cores l|Paczvhski et al.ll2007l ). 

A possibility of existing cool contact binaries with 
secondaries in a more advanced evolutionary stage than 
primaries has been mentioned in a number of papers 
jTapia fc WhelanI Il975l : ISarna fc Fedoroval Il989l : lEggletonl 
Il996l ) but such systems were considered to be exceptions 
from the general rule stating that W UMa stars have not re- 
versed mass ratios. An evolutionary model with mass ratio 
reversal being a necessary co ndition f o r forming a W UMa 
type binary was proposed bv IStepiehl l|2004 l2006al lbh. The 
starting configuration for such a model is a close detached 
binary with a initial orbital period of a couple of days. The 
binary loses AM via a magnetized wind from both com- 
ponents. The AM loss rate varies approximately as M~^, 
where M is a stellar mass (jCazeas fc Stepiehlliooil ). This 
dependence is very similar to the mass dependence of the 
MS evolutionary time scale. Because of this coincidence the 
primary is expected to be close to, or just beyond termi- 
nal age MS (TAMS) at the time when the Roche lobe de- 
scends onto its surface and the mass transfer to the sec- 
ondary begins. Conservative, or nearly conservative mass 
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transfer makes the orbit shrink, which accelerates the trans- 
port rate. The process stops after majority of the primary's 
mass is transferred and the mass ratio reversed (i.e. the for- 
mer secondary component becomes now the primary). De- 
pending on the detailed values of the binary parameters and 
the mass transfer process, a con tact binary or a very short 
period Algol-type star emerges l|Stepiehll2006ah . In the lat- 
ter case an additional AM loss via the magnetized wind is 
needed to turn the Algol into a contact binary. After receiv- 
ing hydrogen-rich material from the secondary, the primary 
moves upward toward ZAMS (corresponding to its larger 
mass), whereas the secondary becomes substantially over- 
sized for its (lower) mass. Each star separately is in thermal 
equilibrium while filling its critical Roche lobe. Further evo- 
lution of the binary is governed by a self-regulating mech- 
anism with two processes acting in the opposite directions: 
evolutionary expansion of the secondary, followed by mass 
transfer to the primary, which leads to widening of the orbit 
and orbital AM loss tightening it. As a result, a contact con- 
figuration is maintained with a slow (on a nuclear evolution 
time scale) net mass transfer from the secondary to primary 
l|Gaze as fc Stcpich 2008) until coal escence occ urs when the 
mass ratio reaches a critical value (|Rasidl 19951 ). 

The model explains the Kuiper paradox, i. e. the exis- 
tence of an equilibrium configuration with component sizes 
required by the Roche geometry. It does not, however, ex- 
plain the observed properties of light curves. The energy 
transport between the components is not a part of that 
model. It is an independent process which will be consid- 
ered in detail in the present paper. 

The paper is organized as follows: basic assumptions 
and equations governing the mass and energy flow are given 
in Sect. 2. The stream structure across and along the flow is 
discussed in Sect. 3, together with application of the derived 
formulas to two actual contact systems. Section 4 describes 
the global reaction of both components to the circulation 
and Sect. 5 discusses the results of the paper. It also sum- 
marizes the main results. 



2 EQUATIONS AND ASSUMPTIONS 

The basic Eulerian equations of fluid flow in a frame of ref- 
erence rotating with the binary are: 
continuity equation 



0, 



momentum equation 

— -I- (v • V V = - Vp -I- i^V V 
ot p 

and entropy equation 



Vip — 217 X V , 



| + (v.V). 



= -V • F + dviB, 



(1) 



(2) 



(3) 



Here p,p, s and v denote gas density, pressure, entropy 
and velocity, ip is the gravitational (plus centrifugal) poten- 
tial, v is kinematic viscosity coefficient, dvisc is the viscous 
term, F is the radiative energy fiux and is angular velocity. 
In the momentum equation the viscous term with V(V • v) 
has been left out as unimportant (see below). 
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Figure 2. The reference frame and geometry of the flow in the 
meridional plane of the secondary component. The stream flows 
away, perpendicularly to the plane of the paper, in the (/>-direction. 
Scales for the transition layer depth and the flow thickness have 
been exaggerated. 



Assuming that the life-time of the cool contact configu- 
ration is of the order of one or more Gyr, we consider a sta- 
tionary solution of the above set of equations on a time scale 
longer than the thermal time scale. We put, then, d/dt = 0. 
We assume that a steady mass fiow from the hotter primary 
to the secondary exists in the common envelope. The as- 
sumption of steady motion eliminates any dynamical insta- 
bilities, generation of acoustic waves or any other processes 
taking place on a very short time scale. After leaving the 
neck between the components the mass stream quickly as- 
sumes the azimuthal motion with hydrostatic equilibrium 
condition in the vertical and meridional directions satisfied. 

Let us first consider an inviscid motion. The inertial 
term on the left hand side of Eq. (2) is of the order of 
v'^/l = /2nRaec, as the flow driving force results from 
the azimuthal pressure gradient. It can be compared to the 
Coriolis term of the order of 2^v. Their ratio is called the 
Rossby number Ro 



Ro = 



27ri?,, 



-/2nv : 



(4) 



where Vcq is the equatorial velocity of a rotating star. 
In a typical W UMa type star Vcq ~ 100 — 150 kms^^, so as 
long as the fiow is not highly supersonic the Rossby number 
is small, i.e. Ro « 1, and we can safely neglect the inertial 
term. The equation of motion becomes then 

-'Vp = 2n xv + Vip, (5) 
P 

where v has only one, non-vanishing azimuthal com- 
ponent V. The above equation describes the so called 
geostrophic fiow in which the meridional component of the 
Coriolis force balances the lateral pressure gradient, per- 
pendicular to the direction of motion. Such flows are often 
observed in the terrestrial atmosphere when air circulates 
around a high or low pressure center. In our case, the Cori- 
olis force conflnes the mass flow to the equatorial region of 
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width depending primarily on the flow velocity (see Sect. 3.4 
below). 



3 THE STREAM STRUCTURE 
3.1 Lateral equilibrium 

We introduce a spherical system of coordinates with the 
origin at the secondary star center (r, (p) and we neglect 
the deviation of equipotential surface shapes from spheres 
(Fig. 2). Because the stream is symmetric relative to the 
equator, we consider only a hemisphere 55 ^ 7r/2. Far 
from the neck matter flows along the equator of the sec- 
ondary in a belt with a half- width Aa; = i?sec(7r/2 — i?f(r)), 
where i?f (r) is the depth dependent polar angle of the flow 
boundary. In general, pressure inside the flow is a function 
of all three coordinates but within our assumptions it is a 
slowly varying function of the azimuthal direction so it can 
be assumed constant when considering equilibrium in the 
radial and meridional direction. The equilibrium condition 
results from Eq. (5): 
meridional component 

= 2pkCsncos^, (6) 

r ov 



and the radial component 

dp dip 
— — 2pkcs^l sin ■& + p—— . 
or dr 



(7) 



Here Cs and k are the sound velocity and the Mach num- 
ber of the flow, respectively. We assume fc ^ 1 throughout 
the paper. Eq. (6) can be integrated 



p(i9) = po — 2rpkcsQ,{l — sin i9) . 



(8) 



where po is the pressure at the equator. The ambient 
pressure under the secondary's photosphere is equal to Pscc, 
so the boundary will be reached when p(i?f) = Psec- This 
gives a condition for i^f 



smff 



Po Psec 



(9) 



2pkCsVcq 

For parameter values characteristic of W UMa stars, 
typical stream widths are around 15° - 60° (see Sect. 3.4). 

Above the photosphere pscc = 0. With a total thickness 
Ar << Rscc the stream layer is thin, so we can replace 
coordinate r appearing explicitly in Eq. (8) by Race- 

If we now include the effects of viscosity, this will pro- 
duce a drift of the stream matter into high astrographic 
latitudes. Molecular viscosity is negligibly small but the con- 
vective viscosity may play a role. We can estimate a possible 
drift resulting from that viscosity 

— "T^conv^conv , (-^^) 



where iconv and Vconv are the mixing length and convec- 
tive velocity near the bottom of the stream. 

Assuming that Zconv is of the order of pressure scale 
height near the bottom of the flow one obtains from the 



model of the convective zone H„ ~ 10 cm and Vc. 



2 X 



7 X lO^'^ in cgs units (|Baker fc Temesvarvl 

1 19661 : ISpruitlll974h . The resulting viscous term in the mo- 
mentum equation is of the order of lyv/Ax^. For v « 0.3cs 
and Ax « i?scc/2 this term is equal to 2 x 10~^ cms~^. 



It takes about 3 d for a given flow element to encircle the 
secondary assuming a typical flow velocity about 10 times 
slower than v^q, see below. The total drift distance of the 
stream in the meridional direction is equal then to ~ 2 x 10'^ 
km which is negligible compared to the stream width. We 
conclude that viscosity can be neglected altogether when 
considering dynamics of the mass flow around the secondary. 

3.2 Vertical equilibrium 

A typical W UMa type binary consists of two unequal 
mass components possessing convective envelopes. Except 
for stars with spectral type around FO where the convection 
zone is very thin and mostly super-adiabatic, its dominant 
part is stratifled nearly adiabatically i. e. according to the 
pressure-temperature relation 



(11) 



where T is temperature, a = 2.5 for fully ionized gas and 
K = const, throughout the considered part of the convec- 
tion zone. The parameter K is called the adiabatic constant 
and its value depends on the speciflc entropy of the matter 
in the convective zone. Realistic models of the convective 
envelope take into account all important effects inffuenc- 
ing its stratification, like variability with depth of A", a, and 
speciflc heats as well as d e viations from strict adiabaticity 
iBaker fc Temesvarvl [l966l : ISpruitlll974l '). However, most of 
these effects take place in the uppermost layers of the con- 
vection zone where the pressure is low. Because we are inter- 
ested in the pressure equilibrium at some depth below the 
stellar surface, a detailed structure of these shallow layers 
is unimportant. As the accurate models of the convection 
zone of solar type stars indicate, the stratification becomes 
very close to adiabatic already for temperatures higher than 
about 10* K. We adopt the adiabatic stratification for the 
convective envelopes of both components with different val- 
ues of the adiabatic constant. 

Matter flowing from the primary is characterized by the 
same adiabatic constant Tfpr as the convection layer of the 
primary. The convection layer of the secondary is described 
by ifscc- In general, A"scc < -?fpr. Figure 3 gives an example 
of two pressure-temperature relations corresponding to two 
adiabatic constants differing by 0.5dex. The upper line on 
each diagram gives the p — T relation for the stream and the 
lower line for the ambient convection layer of the secondary. 

Vertical hydrostatic equilibrium condition requires con- 
tinuity of the pressure at the boundary between the stream 
and the underlying convective layer not affected by the flow. 
Because the secondary convection zone has specific entropy 
lower than the stream, a discontinuity in entropy and tem- 
perature must occur at this boundary (a vertical segment of 
the solid line in Fig. 3 top) . A steady state model with a con- 
tact d iscont inuity of this kind was considered by Shu et al] 
(j 19761 . 1 19791 ) . They assumed that matter streaming from the 
primary fully covers the secondary component of a contact 
binary. They tried to build a model in which the nuclear 
energy of the secondary can fiow outwards in spite of the 
discontinuitv. iHazlehurst fc Refsdal (1978. ) argued, however, 
that such a model contradicts the second law of thermody- 
namics because the nuclear energy of the secondary cannot 
flow from the low to the high entropy medium above. In- 
stead, the nuclear energy will heat the convection zone on 
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Figure 3. The adiabatic p — T relations for the stream matter 
(upper lines) and the secondary's convection zone (lower lines). 
The condition of the continuous pressure across the contact sur- 
face between the stream and the unaffected convective layer be- 
neath requires a temperature discontinuity (a solid line, top) . The 
formation of a radiative transition layer with a flat temperature 
gradient between the stream and the deeper unaffected layers re- 
sults in a continuous transition both in pressure and temperature 
(a solid line, bottom). 



a thermal time scale rising its specific entropy to the value 
characteristic of the blanketing matter, just as assumed in 
the TRO model. 

The situation considered in the present paper is, how- 
ever, entirely different. The streaming matter covers now 
only a part of the secondary. There is no need to introduce 
a contact discontinuity. The hot stream matter heats the 
layer lying underneath and lowers the temperature gradient 
there. As a result, the convective energy flow from below is 
inhibited. In steady state a radiative layer is formed with 
a low temperature gradient and specific entropy increasing 
outwards from the value characteristic of the secondary con- 
vection zone up to the value characteristic of the stream gas. 
We will call it a transition layer. Such a layer with the tem- 
perature gradient equal to zero is shown in Fig. 3 (bottom) 
as a horizontal segment of the solid line. In real stars the gra- 
dient will probably be not equal exactly to zero so the tran- 
sition layer will be thicker than shown. The heat transport 
is completely, or nearly completely, blocked in the transition 
layer which acts as an insulator. The energy blocked by the 
stream will be redistributed over the whole convection zone 
of the secondary and ultimately radiated away from the po- 
lar regions not covered by the stream matter (see Section 
4). 

To consider the vertical structure of the stream in more 
detail we need to calculate po — Psec with depth. Because 
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Figure 4. Tiic vertical pressure distribution across the stream to- 
gether with the transition layer (broken line), in the unperturbed 
secondary convection zone (dotted line) and their difference (solid 
line), balanced by the meridional component of the Coriolis force 
(see text). The bottom of the transition layer is at z = and 
the bottom of the stream (expected to coincide with the Roche 
critical surface) is at z 4 x 10^ km. 



the depth of the layers in question is much lower than the 
stellar radius we can use the plane parallel approximation 
with the coordinate z replacing r (Fig. 2). Let us assume 
that the secondary convection zone remains unaffected by 
the stream below the bottom of the transition layer and 
we put the reference level z = there . The pressure is 
constant over the whole equipotential surface lying at that 
level. We denote it by pb- Within the isothermal transition 
layer pressure varies as 



p{z) = Pb exp{-z/Hp), Hp 



Q^z^Zo, (12) 



where /i is the mean molecular weight, H - mass of the 
hydrogen atom, ks - the Bolzmanu constant, To - temper- 
ature of the layer, Qes - the effective gravity and Zo corre- 
sponds to the bottom of the stream, which agrees with the 
top of the transition layer. The pressure is continuous at the 
interface between the transition layer and the adiabatically 
stratified stream where temperature varies as 



T{z)=To-^{z~Z„), Z^Zo: 

Cv 



(13) 



where Cp is specific heat at constant pressure. When 
calculating the T[z) relation wc neglect variability of the ef- 
fective gravity. With T{z) known, the vertical pressure strat- 
ification is obtained from Eq. (11) where K = KgtT- As it 
will be shown below, the entropy of the stream matter varies 
very little during the flow around the secondary, so we can 
put ifstr ~ everywhere. 

It is reasonable to assume that in steady state the tran- 
sition layer will be partially dragged by the stream in such a 
manner that a gas layer lying immediately below the stream 
moves with a stream velocity and deeper layers move slower 
until zero velocity is attained at the bottom of the tran- 
sition layer. We arbitrarily assume that velocity inside the 
transition layer varies as 



v{z) = kc, 



Pb - p{z) 

Pb - Pt 



(14) 
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where pb and pt are the densities at the bottom and at 
the top of the transition layer. Thus defined velocity vanishes 
at the bottom and reaches the stream velocity at the top. 

Far from the stream the convection zone of the sec- 
ondary remains unperturbed with the adiabatic stratifica- 
tion. Neglecting a difference between Cp in the stream and 
in the ambient convection zone we have 



T{z) = To 



5cff 



(15) 



Pressure stratification is obtained again from Eq. (11), 
except that now K — Kscc- 

Figure 4 shows an example of the vertical pressure dis- 
tribution inside the transition layer and the stream with 
depth equal to 1 % of the primary's radius (broken line), 
in the unperturbed secondary convection zone (dotted line) 
and their difference (solid line). As it is seen from the fig- 
ure, the stream extends above the surface of the secondary 
and produces an equatorial bulge held by the meridional 
component of the Coriolis force (see also Fig. 5). 



3.3 The thermal stream structure along the flow 

With the depth and width of the stream specified, the mass 



From Eqs. (17) and (19), the limiting case of large heat 
capacity, dQ » dL, corresponds to the condition 



CpMTAr » 



2TvaTp^Rscc 
kcs 



(20) 



If the radiated energy is small compared to the total 
heat capacity of the stream, its effective temperature will 
change very little from its initial value Tpr so we can put 
Tstr ~ Tpr- In case of low heat capacity the effective temper- 
ature of the stream will approach the effective temperature 
of the secondary and the stream matter will mix together 
with its convection zone. 

Assuming that the energy outflow from the stream in- 
fluences uniformly its entropy, i. e. that the stream specific 
entropy is constant with depth, we can calculate its decrease 
from Eq. (3). In case of the large heat capacity we obtain 



As : 



-i-irRsccCrTp^ 

kCsMTAr 



(21) 



In case of the low heat content, Tpr in Eq. (21) should 
be replaced by the properly defined average stream effective 
temperature. 

The entropy decrease can be compared with entropy of 
the primary convection zone 



Fp = 2Ax kcspdr ^ 2AxkcsMAr , 
'o 



(16) 



where Mav is the column mass above the level Ar. The 
stream also carries thermal energy which is partly radiated 
away during its flow around the secondary. We will now 
discuss the variation of the stream thermal structure along 
the flow. 

Let us consider a slice of matter perpendicular to the 
flow with thickness i?sccd<j!>. Its total heat (internal energy 
plus enthalpy) capacity dQ is equal to 



dQ = 2AxRs. 



where 



MTat 



pTdr. 



pCpTdr = 2AxR,ccCpMTArd(j}, (17) 



(18) 



The matter flowing around the secondary radiates en- 
ergy at a rate aT^tr per unit time and area, where Tstr is 
the effective temperature of the stream. In general, Tatr will 
decrease along the flow from its initial value, assumed to 
be close to the effective temperature of the primary, Tpr, 
and its precise variation could be determined by solving the 
transfer equation within the stream. To avoid cumbersome 
calculations we discuss two limiting cases: first, when the to- 
tal heat capacity of the stream is much higher than the total 
energy radiated away, and then, when the heat capacity is 
comparable to the radiated energy. 

The slice has a radiating area 2Aa;i?sccd0 and emits 
dFrad ~ 2(7 Tstr Ax_Rsccd(^ energy per unit time. It takes a 
time equal to 2nRscc/kCs to encircle the secondary. The total 
energy radiated away during this time is 



dL = 



47ro-r,t,.Axi?; 



kcs 



(19) 



(22) 



After return to the primary, the stream matter sinks in 
its convection zone. 



3.4 Numerical example: AB And 

To obtain a numerical estimate of the parameters describ- 
ing an actual stream in a cool contact binary we apply the 
above derived equations to a typical W-type contact bi- 
nary. As an example, we take AB And. It has the orbital 



period Porb = 0.33 d, component masses 



1.04M, 





and Msec = 0.6 0Mq, and volume radii 7?pr = l.O2_R0 and 
-Rsoc = O.78i?0 l|Baran et al.ll2004l ). These data give Uoq = 
120 kms~^ for the secondary compone nt. The star has a 
relatively low fill-out factor / = 0.05 (|Baran et al.ll2004l ). 
which corresponds to about 1 % of the primary star radius. 
Assuming that the equal pressure equipotential surface is 
close to the inner critical Roche lobe we obtain Ar ~ 0.01 
for the depth of the stream. The sound velocity is equal to 
30 kms~^ at this depth. The width of the stream can be 
calculated from Eq. (9). Figure 5 shows the results. Here 
the angular half-width is plotted as an abscissa for differ- 
ent heights above the bottom of the transition layer. The 
broken heavy line on the right indicates the approximate 
photospheric level of the secondary. The diagrams give the 
stream width for three values of Ar/Rpi equal to 0.005 (top), 
0.01 (middle) and 0.03 (bottom) and three values of k, as 
indicated. The width was calculated up to the level where 
temperature reached 10 000 K because for lower temper- 
atures the adiabatic model breaks down. As we see, the 
largest width of the stream is close to the photosphere of 
the secondary but the Coriolis force produces a bulge above 
the stellar equator with a height of the order of the stream 
depth. As expected, a lower stream velocity results in its 
increased width and in the limit of zero velocity the stream 
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Figure 5. The width of tiie stream as a function of lieigtit above 
the bottom of the transition layer for three values of k describing 
the ratio of stream to sound velocity and for three depths of 
the stream: 0.5 % (top), 1 % (middle) and 3 % (bottom) of the 
primary's radius. A broken heavy line on the right indicates an 
approximate level of the secondary's photosphere. Note a weaJi 
dependence of the stream width on its total depth. 



would cover the whole secondary. For flow velocities com- 
parable to the sound velocity, i. e. for k between 1 and 0.3, 
the stream hcis a half- width between 15° and 30° (Fig. 3), 
and covers between 1/4 and 1/2 of the stellar surface, re- 
spectively. 

A schematic side view of AB And with a stream covering 
50 % of the secondary's surface is shown in Fig. 6. 

The stream with a depth of 0.01 carries enough mass 
and energy to be in the regime of large heat capacity. To 
see it we substitute the following values to Eqs. (17) and 
(20): Ar = 7 X 10^° cm, Cp = 10^ (in cgs units), Tpr = 
5590 K (see next section) and k = 0.3. The column mass 
and the integral MTj^r are calculated from the convec- 
tion zone model. We obtain Mat- = 3.5 x 10^ gcm~^ and 
MT Ar = 1.7 X 10^° gKcm"^. The resulting mass transfer 
rate is Fp = 5x 10"" M© /year. The left hand side of Eq. (19) 
gives 1.7 X 10^^ (in cgs units) and the right hand side gives 
4 X 10^® in the same units. The resulting ratio of both num- 
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Figure 6. A side view of AB And with a stream covering a half of 
its surface. Broken lines indicate inner and outer critical surfaces 
and the solid line shows the binary surface. 



bers is equal to 2.4 x 10~ which indicates that the non- 
equality given by Eq. (19) is very well fulfilled. The numbers 
become comparable only for stream depths less than 0.5 % of 
the stellar radius. Figure 5 (top) shows the calculated width 
of the stream with the depth of 0.5 %. It is somewhat sur- 
prising that its width is not much different from the widths 
of much more massive streams shown in Fig. 5 (middle) and 
(bottom). It is clear that only extremely marginally contact 
binaries with overfill factors less than a couple of hundredths 
will have streams with surface temperature approaching the 
temperature of the uncovered part of the secondary compo- 
nent. 

The change of specific entropy of the stream in AB And 
can be calculated from Eq. (20). With the adopted param- 
eters we obtain As = —1.2 x 10^ ergg~^K~^. According to 
Eq. (21) the initial value of the specific entropy is equal to 
4.59x10* in the same units so the decrease of entropy during 
the flow around the secondary is small indeed. It is possible 
then that the stream returning to the primary component 
will not sink immediately after passing the neck but it may 
move over a substantial fraction of the primary's circumfer- 
ence before it plunges in the convection zone. 

To summarize, we see that the Coriolis force con- 
fines the near sonic stream to the equatorial belt of width 
±(15° — 30°) around the secondary. The stream covers 25- 
50 % of the stellar surface. For a primary overflowing its 
Roche lobe by --^ 1% or more the stream radiates a very 
small fraction of its total heat capacity which results in 
its nearly constant surface temperature equal to the pri- 
mary's temperature. Substantially shallower streams, possi- 
ble in contact binaries with an extremely marginal contact, 
will show a temperature decrease from the initial value close 
to the primary's temperature on the trailing hemisphere of 
the secondary, down to the value close to the secondary's 
temperature on the leading hemisphere of the star. 

After return to the primary, the stream will sink in 
its convection zone due to the lower entropy. The Corio- 
lis force deflects the returning stream to the opposite side of 



8 K. St§pien 



0.4 r 



0.3 7 

p. 

X 0.2 7 

<i : 




Figure 7. A fraction of the primary component luminosity radi- 
ated away by the stream flowing around the secondary component 
as a function of mass ratio. The curves correspond to three values 
of coverage of the secondary by the stream: full coverage (dotted 
line) , half-coverage (solid line) and quarter-coverage (broken line) . 

the neck, compared to the stream flowing towards the sec- 
ondary. The deflection should reduce an interaction of both 
flows although it will not separate them completely. The re- 
turning flow will collide with matter stre aming towards the 
neck (see Figs. 1 and 2 in lOka et aLll2002l ). However, the col- 
lision will not be head-on, as obtained bv lMartin fc DavevI 
||1995[1 who simulated a two dimensional motion in an equa- 
torial plane of a contact binary. Geometry applied by the 
authors prevents matter from crossing the plane containing 
centers of both stars and the rotation axis, so the stream 
cannot move to one sid e of the neck, as shown by three di- 
mensional calculations (|Oka et al.ll2002l ). 



4 GLOBAL REACTION OF THE BINARY TO 
THE MASS AND ENERGY FLOW 

4.1 The influence of the stream 

As it was mentioned earlier, the stream acts like a hot blan- 
ket covering part of the secondary and blocking the energy 
flow from below by rising the temperature of the layer im- 
mediately beneath. As a result, a transition radiative layer 
is formed with a decreased temperature gradient in which 
the specific entropy increases from the bottom value char- 
acteristic of the secondary convection zone to the top value 
characteristic of the stream. 

What is the global reaction of a star when a fraction 
of its surface cannot radiate away the stellar core energy 
due to an obstacle blocking the energy outflow? The phys- 
ical nature of the obstacle is not so important as long as 
it occurs in the outermost stellar layers. One example of 
such an obstacle partly blocking the energy flow from be- 
low is a dark starspot appearing on the surface of an active 
star. Another example can be seen in close binaries when 
a hot companion irradiates a part of the surface of a cool 
component. Both e ffects have been discussed at length in 
the literature (e. g. iRitter et al]|2000l , [Spruit fc Weissill986l 
and references therein). The blocked energy flux is redis- 
tributed inside the convection zone and re-radiated by the 



unperturbed part of the stellar surface. The detailed mod- 
els show that the specific entropy of the convection zone is 
somewhat increased and so is the stellar radius. For stars 
with deep convection zones, or fully convective, the nuclear 
energy production may also be affected. Assuming a steady 
state in which a constant fraction of the stellar surface is cov- 
ered by spots over a time lo ng compared to the the rmal time 
scale of the convection zone, ISpruit fc Weiss ('1986') obtained 
the following expression for the equilibrium stellar radius 

Re{fps)=Ro{l~fpsf, (23) 

where Re and Ro denote the equilibrium radius of the 
spotted and unspotted star, fps is the fraction of the stel- 
lar surface permanently covered by spots an d d ~ —0.1. 
An ess entially identical expression was found bv lRitter et al.l 
(2000') for stars irradiated by hot companions. For fps = 0.25 
and 0.5 the equation gives R^ = 1.03-Ro and 1.07-Ro, respec- 
tively. In stars with not too deep convection zones the nu- 
clear energy production remains unaffected. The reduction 
of the effective radiating surface is exactly compensated by 
an increase in surface temperature of the undisturbed part 
of the star 

rcff = (l-,fp.)"^Tcff, (24) 

where T^g and Tefj are the equilibrium and undis- 
turbed effective temperature. More accurate values of the 
equilibrium radius an d temperature depend on stellar mass 
ISpruit fc Weissll 19861 ). 

The primary component suffers from additional cool- 
ing of its surface by the stream encircling the secondary. 
The stream feeds the primary with the lower entropy mat- 
ter which is mixed with the rest of its convection layer. As 
a result, the specific entropy of the convection zone is some- 
what lowered and so is the stellar radius but the expected 
changes should be close to those occurring in the secondary 
component, i.e. they should not exceed one, or at most a 
few percent. The primary, slightly undersized for its mass, 
will appear closer to ZAMS in the HR diagram. Neglect- 
ing its radius change we can obtain the modified effective 
temperature of the primary 

TpT = (1 - ^) , (25) 

where Lpr is the luminosity of the primary component 
and ALpr is a fraction of this luminosity radiated away 
by the stream flowing around the secondary component 
(Fig. 7). 

Now let us calculate the expected values of the thermal 
parameters of AB And. The effective temperature of a single 
one solar mass star is about 5800 K. That should be close to 
the temperature of the undisturbed primary. Assuming that 
secondaries of W UMa type stars arc evolved stars with hy- 
drogen depleted in their stellar cores (Stcpieri 2004, 2006j|) 
we obtain the following values of the global parameters for 
the unperturbed secondary of AB And: T^g = 4570 K, and 
L/Lq — 0.24. The data were taken from the evolutionary 
models of stars with masses equal to the present masses of 
the secondaries. This is not fully consistent with the present 
considerations because, according to the accepted model, the 
secondaries were originally more massive and their present 
masses result from mass exchange. It is hoped, however, that 
differences between the adopted and fully consistent models 
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are of secondary importance. Let us assume, as an example, 
that the stream covers 50% of the secondary and that its nu- 
clear energy production rate is unaffected by this. According 
to Eq. (23) the steady state effective temperature of the sur- 
face fraction not covered by the stream rises to 5340 K. The 
primary transfers about 16 % of its luminosity (see Fig. 7) to 
the secondary so its actual effective temperature, resulting 
from Eq. (25), drops to 5590 K. The effective temperature 
of the stream is of course the same. The "final" surface aver- 
aged temperature of the secondary is eq ual to about 5480 K , 
which is nearly identical as obtained bv lBaran et al.l (|2004| ) 
from observations who give Ts^c = 5500 K. Such a small dif- 
ference between the predicted and observed temperature is 
very likely fortuitous, taking into account a large observa- 
tional uncertainty of the observed value (see below). How- 
ever, the observed primary's temperature, equal to 5140 K, 
is substantially lower than resulting from the present model. 
A possible explanation of this discrepancy is discussed in 
the next Section. 



4.2 The W-type phenomenon and photospheric 
starspots 

The observed temperatures of components of a W UMa type 
binary be long to the least accurate ly determined stellar pa- 
rameters l| Yakut fc Eggletonll2005h . The primary's temper- 
ature is usually determined from the spectral type which is 
uncertain due to very broad spectral lines, or from the pho- 
tometric index. It is kept constant during the light curve 
modeling. The secondary's temperature results from model- 
ing of the light curve, in particular from the relative depth 
of both minima, which gives a ratio of both temperatures. 
Such an analysis gives only a rough estimate of the sur- 
face averaged temperatures and is insensitive to their possi- 
ble variations across the stellar surface. Regarding the rel- 
ative temperatures of the components, cool contact bina- 
ries are divided into two types: W-type when secondaries 
look hotter than primaries an d A-type when the opposite 
takes place (iBinnendiikI Il970h . A typical temperature dif- 
ference between the components of genui ne contact binaries 
is of the order of a few hundred kelvin l|Yakut fc EggletonI 
I2OO5I ). A few stars alternating between W-type and A-type 
are known with a temperature difference always staying close 
to zero. The present model gives AT = Tpr — Tscc equal to 
110 K for AB And. Several A-type stars are known with 
AT close to this value. However, AB And is of W-type 
and this phenomenon cannot be explained within the en- 
ergy transfer mechanism alone. No such mechanism, obey- 
ing basic laws of physics, can transfer energy from cooler to 
hotter medium in an isolated system. So, the most popu- 
lar explanation of the W-type phenomenon assumes that 
dark, cool spots cover a substantial fraction of the pri- 
mary's surface l|Ruc inski 1993.) . The temperature of its spot- 
free part may be higher than the secondary's temperature 
but the apparent surface averaged temperature is lower. W 
UMa-type stars have a very high level of chromospheric- 
coronal a ctivity fCruddace & Duprcc "l98^; iRucins k1 1 19851 : 
IVilhu fc Walter 1987: Stepieh ct al. 2001. ') and variations of 
the light curve, observed on time scales of decades, indi- 
cate that they also possess da rk starspots covering a variable 
fract i on of the stellar surface llGazeas. Niarchos fc Gradoulal 
l2006l : ICoughlin. Dale fc WilhamonlbOOSl ). Perhaps the most 



extreme example of light curve vari ability due to a vari- 
able c overage by spots can be seen in lRucinski fc Paczvnskil 
(l2002h . 

We can estimate how large should be the fraction of the 
primary's surface covered by dark spots to lower its average 
temperature to a required value. Assuming, for simplicity, 
that umbrae have surface brightness of 1/4 of the undis- 
turbed photosphere (it is so in the solar case) and that they 
dominate the spots we have 



r^, = (1-0.75/. )^/*rp 



(26) 



where Tp^ is the "final" primary's temperature modified 
by spots and fs is the spot coverage. It is assumed here that 
spots appear and disappear on a time scale short compared 
to the thermal time scale of the convection zone. When the 
spots appear, they simply block the respective energy fiux 
which does not reappear elsewhere. In other words, we con- 
sider now only the variable component of the spot coverage. 
If there also exists a permanent spot component covering 
a fraction fps of the primary's surface it will infiuence the 
effective temperature as discussed in Sect. 4.1. 

For fs = 0.25 we obtain T^, = 5310 K, hence AT = 
-170 K. When is increased to 0.5, we obtain Tp, = 4970 
K, hence AT = -510 K. This should be compared with the 
observed value AT = -340 K. As we see, the required spot 
coverage is quite high but reasonable. There are a number of 
single, rapidly rotating stars w ith a variab le spot coverage 
within this range, e . g. AB Dor l|lnnis et al.|[2 008'l or BO Mic 
(Wo lter et aLllioOSh . Heavily spotted are also cool detached 
binaries with period s shorter th an one day, like XY UMa 
iPribulla et al.ll200ll ) or BI Get (|Cutispoto et al.ll2003l ). 

The absolute magnitude of AB And, obtained from the 
visual magnitude and the distance is equal to 4.0 5 mag 
jR ucinski fc Duerbeck|[l997[ ). With BC = -0.08 mag (|Allenl 
1197^ ) the binary has the luminosity L ~ I.6L0. Adopting 
L « I.I6T0 for the MS primary of AB And and L ^ 0.251/© 
(see above) for the evolved secondary we obtain IALq for 
the total luminosity of the binary. This is less than observed 
but the difference is not large and can be due to observa- 
tional uncertainties. Note that if the secondary were a nor- 
mal MS star, its core luminosity would be close to O.ILq 
which increases the discrepancy with observations. 



4.3 The case of AW UMa 

Up to recently the spectroscopic observations of contact bi- 
naries were not accurate enough to resolve an additional 
fiow on a rapidly rotating component. High-precision obser- 
vations obtained in the last ye ars by S. M. Rucinski and 
his group changed the situation. iPribuUa fc Rucinskil (|2008| ) 
analyzed very accurate profiles of spectral lines of AW UMa 
- an A-type contact binary with an extreme mass ratio. The 
observations indicate the existence of an equatorial bulge on 
both components, moving with velocity 20-30 km/s relative 
to the rotating frame of reference. The present model pre- 
dicts such a bulge on the secondary star. Although dynamics 
of the returning stream was not considered in the present 
paper it was argued that the returning flow may travel a 
significant fraction of the circumference of the primary be- 
fore it disappears beneath the stellar surface. The returning 
fiow, together with the stream formed on the leading hemi- 
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sphere of the primary, wiU produce a spectral feature similar 
to signature of an equatorial bulge. 

The observations bv lPribulla fc Rucins ki' (2008') showed 
in addition that parts of the line profiles formed away from 
the equatorial regions of AW UMa are too narrow for the 
rigidly rotating components filling the Roche lobes. The au- 
thors suggest that the components may be undersized by 
about 15 percent. However, another explanation is als o pos- 
sible. The flow pattern obtained bv lOka et al.l (|2002h from 
hydrodynamic simulations of a semi-detached binary shows 
the presence of giant eddies surrounding the poles of the 
donor and rotating contrary to the orbital motion (i. e. con- 
trary to the direction of the stream). Such eddies reduce 
rotational broadening of lines formed in the polar regions. 
In a contact binary similar eddies may also occur on the 
secondary. If so, the observed profiles mimic rigid rotation 
of undersized stars. 

IPribuUa fc Rucinskil (|2008| ) stress, however, that such 
profiles are not common and are not observed e. g. in an- 
other contact binary, V556 Oph, in which the spectroscopic 
observations are in agreement with a conventional contact 
model. The difference between line profiles observed in AW 
UMa and V556 Oph may result from a different evolution- 
ary status of both stars and different size of the stream. 
The secondary component of AW UMa is very likely to be 
a highly evolved s tar with a massive heli um core and a very 
tenuous envelope (|Paczvnski et al.ll2007l ). The evolutionary 
computations suggest that the size of the secondary may be 
close to its maximum, or even past it, so the secondary may 
presently undergo evolutionary shrinking. In such a case, the 
stream can be more massive and move faster than in other 
binaries in which secondaries still expand due to evolution- 
ary effects. A contact configuration of AW UMa may still 
be sustained by a slow (net) mass transfer from the primary 
to the secondary. This transfer results in a shortening of or- 
bital period i. e. it acts in the same direction as AML. The 
situation when two operating mechanisms tighten the orbit 
cannot last for long. If so, AW UMa is in an exceptional 
evolutionary state just prior to final merging of both com- 
ponents. It would be very interesting to identify other stars 
with properties similar to AW UMa. Secondaries with little 
or no helium core expand as they evolve. A self-regulatin g 
mechanism works in such binaries (IGazeas fc Stepieri|[200i ). 
The expanding secondary transfers mass to the primary on 
an evolutionary time scale. Mass transfer from a less massive 
to a more massive component widens the orbit. At the same 
time, AML tightens the orbit on approximately the same 
time scale. The whole process contains a negative feed-back: 
too fast secondary expansion results in an increased mass 
transfer, which leads in turn to an appropriate widening of 
the orbit (in spite of AML) and the Roche lobe increase. The 
Roche lobe of the secondary approaches its surface and the 
mass transfer is cut down. If, for some reason, the expansion 
rate of the secondary becomes low, the orbit tightens due to 
AML, the Roche lobe shrinks and the higher mass trans- 
fer rate is restored. Similarly, if AML significantly increases 
(decreases), the Roche lobe varies correspondingly so that 
mass transfer rate increases (decreases) , resulting in an orbit 
which is almost insensitive to such fluctuations. V556 Oph 
and many other W UMa type stars are very likely in such 
an evolutionary state. Their streams will be less violent and 
not so contrasted as in case of AW UMa. 



5 SUMMARY AND DISCUSSION 

5.1 The essentials of the circulation model 

It was shown by IStepienI l|2004 l2006al lbl) that an evolution- 
ary model of a cool contact binary, in which the primary 
component is a MS star and the secondary component is in 
an advanced evolutionary stage with a hydrogen depleted 
core, fulflUs simultaneously two conditions: both compo- 
nents are in thermal equilibrium and their sizes conform 
to the geometrical requirement of the Roche model. Energy 
exchange between the components was not a part of that 
model. The present paper considers the problem of the en- 
ergy transport. 

A common envelope of a contact binary with unequal 
components cannot achieve hydrostatic equilibrium. Un- 
equal heating from below of the base of the common enve- 
lope by the emergent energy fluxes of the components pro- 
duces temperature difference on each equipotential surface. 
A resulting baroclinic structure in the common envelope 
drives large scale mass motions between the components 
on a dynamical time scale. Because the stream transports 
thermal energy together with mass, stellar luminosities are 
redistributed over the common surface. 

Three dimensional hydrodynamic simulations of the 
mass loss in a semi-detached binary demonstrate that a 
stream leaving the donor star is deflected by the Coriolis 
force in the direction of the orbital motion so that matter 
flows at an angle of about 10° to the line joining star centers 
l|Oka et al.l I2OO2I ') . The speed of the stream reaches sound 
speed already in the vicinity of Li point iLubow fc Shul 
I1975I : lOka et all I2OO2I '). After leaving the Li region the 
stream moves down the potential well. If the companion 
is compact enough, the stream misses its surface, encircles 
it and forms a disc. If a radius of the companion is larger 
than a certain critical value, the stream hits its surface. The 
subsequent stream motion is restricted to two dimensions 
along the stellar surface. A contact conflguration is an ex- 
treme case of a large companion filling the same equipo- 
tential surface as the donor. The stream is then forced to 
move from the beginning along the common equipotential 
surface, to encircle the other component and to return to 
donor. The same Coriolis force which deflects the stream, 
prevents it now from spreading up to the stellar poles. Even 
in the presence of turbulent viscosity the spreading is in- 
significant. The stream flows in an equatorial belt with the 
width determined by the stream velocity. Its part is raised 
above the photosphere and forms an equatorial bulge with 
height comparable to the stream depth. 

Assuming that the bottom of the stream is close to the 
inner critical Roche surface we can calculate the global pa- 
rameters of the flow. Typical overfill factors of the W UMa 
type binaries correspond to thickness of the common enve- 
lope of the order of one o r a few percent of the stellar radius 
(|Yakut fc Eggletonll2005h . For the stream depth of one per- 
cent and velocity close to the sound speed the stream has a 
total width of 30°-60° and carries about 5 x lO"** Mo/year 
with velocity of 10-30 km/s. Its heat capacity is a few orders 
of magnitude higher than the energy radiated away during 
the flow around the secondary. As a result, the effective tem- 
perature of the stream and its specific entropy decrease very 
little when it returns to the primary. The returning flow may 
travel a significant fraction of the stellar circumference be- 
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fore it sinks. Only binaries with an extremely marginal con- 
tact in which the common envelope is thinner than about 0.5 
percent of the stellar radius have streams with heat capacity 
comparable to the energy radiated away. The effective tem- 
perature and specific entropy of such a stream approach the 
ambient temperature and entropy of the secondary's con- 
vection zone. 

The stream covering a fraction of the secondary modi- 
fies its thermal structure. It heats the matter immediately 
beneath reducing the temperature gradient. The reduced 
gradient effectively blocks the energy flow from below. A 
similar situation occurs in stars covered by dark spots or 
irradiated by hot companions. The blocked energy is redis- 
tributed inside the convection zone and radiated away in 
the polar regions not covered by the stream. The effective 
temperature of these regions rises accordingly but the stel- 
lar radius hardly changes. Numerical examples show that 
the surface averaged temperatures of both components are 
close to one another as observed in contact binaries but the 
calculated primary's temperature cannot be lower than the 
secondary's. This is a known result; mass and energy trans- 
fer between the components cannot by itself result in the 
secondary's temperature higher than the primary's. The ex- 
istence of A-type binaries can be explained solely by large 
scale circulations but not W-type. As has been suggested 
many times in the past, large dark spots appearing tem- 
porarily on the primary can explain its low surface averaged 
temperature. Numerical estimates obtained in the present 
paper indicate that a substantial fraction of the star (25-50 
percent) must be covered with spots to lower the tempera- 
ture by several hundred degrees, required by observations. 



5.2 Discussion 

Dynamics of the mass transfer between the components of a 
contact binary, conside red in th e prese nt paper, is similar to 
the one investigated bv lTassoull (Il992l). The driving force of 
the circulations discussed bv lTassoul comes also from 

the baroclinic structure of the common envelope and the flow 
is influenced by the Coriolis force. The author applied, how- 
ever, a different boundary condition at the stellar surface. 
Based on the observations available at that time he adopted 
zero velocity and the strictly uniform temperature on the 
stellar surface of the secondary. Yet more recent observations 
show the e xistence of massive flows in some of the W UMa 
type stars (|Pribulla fc Rucinskill200^ . On the other hand, 
very little is known about the temperature distribution over 
the surfaces of both components although nonuniformities 
in their surface brightness are n otoriously invoked t o explain 
the observed light curv e s (e. g . iGazeas et al.ll2006l '). 

lYakut fc E ggleton' ("20051) called attention to a possi- 
ble role of lateral energy transfer connected with differential 
rotation observed on the Sun and solar-type stars. Differen- 
tial rotation results from a coupling between rotation and 
turbulent convection in an axially symmetric star. AM is 
transported in the meridional direction by Reynolds stresses, 
meridional circulations and viscous diffusion (Brown et a l. 
I2OO8I ) . The mechanism considered in the present paper is dif- 
ferent from that. The differential rotation was not included. 
Nevertheless, it may play a role in forming a stream which 
transfers mass and energy between the components although 



more elaborated model of the meridional AM transport in a 
contact binary is needed. 



W UMa type stars show often period variations 
(|Kreiner et al.ll200lf ). Systematic changes with values up to 
P — ±(I0~^ — I0~^) day/year have been detected in sev- 
eral stars and interpreted as a signature of TRO. However, 
recent observations of several hundred contact binaries mon- 
itored in the framework of the program OGLE showed that 
the distribution of period variations of W UMa stars can 
be approximated by a normal distribution with an average 
equal to zero and d ispersion equal to 2.3 x lO"'^ day/year 
l|Kubiak et al.l l200dl. The observations are of very high ac- 
curacy but they are extended only over 13 seasons. The anal- 
ysis of the observed period variations suggest that they have 
a random character with predominantly low values of P re- 
sulting in time scales substantially longer than the stellar 
thermal time scale. Such a distribution excludes TRO as 
a primary source of the variations. A much more probable 
reason for the observed variations is connected with possible 
fluctuations of the mass transfer between the components. 
The stream flowing from the primary to the secondary and 
back carries I0~'^ — fO~*M0/year. A relative fluctuation of 
this flux at the level of 10~^ — f O"** can produce the observed 
period variations. Fluctuations of that order are expected 
as a result of magnetic activity cycles operating on W UMa 
stars ( ApplcR'atc 1992)- The existence of activity cycles have 
been suggested from the analysis of a long-time photomet- 
ric behavior of heavily spotted stars. The average magnitude 
of those stars varies on time scales of several decades. This 
would be the expected time scale of period variations in con- 
tact binaries. 



The present model of large scale circulations supple- 
ments the evolutionary model of W UMa type stars in which 
contact binaries are past mass exchange with mass ratio re- 
versal ( jStepieri..2004 , 2006a .b) i. e. they are in an evolution- 
ary state similar to low mass Algols. The basic difference 
between W UMa type stars and low mass Algols is that 
the latter stars have more orbital AM and longer periods so 
only the secondary fills its Roche lobe. Binaries with lower 
orbital AM form contact binaries in which the overflow of 
the critical equipotential surface by both components drives 
large scale circulations encircling the whole secondary and 
a large fraction the primary. Matter transported by the cir- 
culations to the secondary (with the mass flux of the order 
of fO^** — fO~*MQ/year) returns back to the primary but 
a fraction of thermal energy associated with that matter is 
radiated away during the flow around the secondary. It was 
argued that the expected flow velocities attain a substantial 
fraction of the sound velocity. Such values were assumed 
when considering details of the flow. To determine an ac- 
curate value of the flow velocity numerical hydrodynamic 
simulations are needed. Particularly important in this re- 
spect is the fate of the returning flow. If it behaves like a 
waterfall, the expected flow velocities can be close to sound 
velocity. If, instead, it collides with the primary's matter so 
that the high pressure wave moves back along the flow, the 
driving force will be reduced and the flow velocity lowered. 



12 K. St§pien 
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